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SIMILrrUDE  MODELING  OF  INTERNAL  GRAVITY  WAVE  SPECTRA 
Edmond  M.  Dew  an 
PL/(OPS),  Hanscom  AFB 


Abstract.  Internal  gravity  waves  are  present  in  both  the  ocean 
and  atmosphere.  The  oceanic  gravity  wave  {jower  spectra 
have  been  modeled  semi-empirically  by  Garrett  and  Munk 
(1972, 1975)  and  shown  to  have  "universal"  properties.  Their 
model  seems  to  be  useful  also  for  atmospheric  waves. 
Recently  a  similitude  type  of  "saturation"  model  for  vertical 
wave  number  horizontal  velocity  fluctuation  spectra  has  met 
with  empirical  support  in  the  atmosphere.  The  present  paper 
employs  a  further  similitude  argument  related  to  "cascade" 
properties  which  would  apply  to  the  frequency  and  horizontal 
wave  number  spectra.  Testable  consequences  are  exhibited. 
If  validated  the  model  would  provide  useful  insight  into  the 
physics  of  these  waves. 

Introduction 

The  power  spectral  densities,  or  PSD's,  of  oceanic  internal 
waves  were  shown  by  Garrett  and  Munk  (1972,  1975)  to  have 
certain  universal  characteristics.  They  also  created  a  ^emi- 
empirical  model  (called  here  GM  72,  75)  for  these  PSD’s. 
Internal  waves  have  also  been  measured  in  the  atmosphere, 
and  Dewan  (1979)  suggested  that  GM  72, 75  may  be  relevant 
to  this  case  as  well.  Van  Zandt  (1982)  has  made  a  convincing 
case  that  this  is  indeed  correct.  Subsequently,  a  striking 
example  of  universality  was  found  by  Dewan  et  al  (1984)  for 
the  vertical  wave  number,  kj,  PSD  for  horizontal  stratospheric 
wind  fluctuations,  u  (Figure  1 )  and  this  demanded  physical 
explanation.  One  based  upon  "saturation"  due  to  turbulent 
instability  was  proposed  by  Dewan  and  Good  (1984),  (1986), 
extended  by  Smith  et  al  (1987)  and  tested  by  numerous 
authors  by  means  of  radar  and  lidar.  Tsuda  (1989)  sum¬ 
marized  the  good  experimental  support  of  this  theory.  Up  to 
that  time,  the  model  described  the  data  very  well.  Most 
recently  Hines  (1991 )  has  proposed  a  new  physical  approach 
which  promises  to  be  very  useful.  It  contradicts  the  previous 
picture  of  "turbulent  saturation",  but  it  does  not  seem  to 
contradict  the  picture  to  be  presented  below. 

The  purpose  of  the  present  paper  is  to  extend  the  similitude 
approach  to  the  remainder  of  the  spectra  of  GM  72,75.  It 
presumes  that,  unlike  the  k^  PSD's,  the  frequency,  to,  and 
horizontal  wave  number,  kjj,  PSD's  are  determined  by  local 
wave  interactions  which  give  rise  to  a  "cascade"  of  energy 
through  time  and  space  scales.  In  a  way,  this  is  like  the  well 
known  turbulence  cascades,  but  certain  wave  properties  are 
retained.  The  final  picture  emerging  is  "wave-like  turbulence" 
of  "turbulence-like  waves". 

Piior  Indications  of  a  Wave  Cascade 

When  Crooks  et  al  (1968)  presented  their  horizontal  wave 
number  PSD’s  of  stratospheric  velocity  fluctuations,  the 
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Fig.  1.  The  Vertical  Wave  Number,  k^,  PSD  of  the  Vector 
Horizontal  Velocity,  u,  in  the  Stratosphere,  V|ty  (k^)  = 

Vu  zonal  (•‘z)  +  Vu  meridional  •’^sed  on  data  in  Dewan  et  al 
(1984).  A  cubic  trend  was  removed.  This  particular  figure 
has  not  been  previously  published.  Note;  These  spectra  are 
not  "normalized".  It  is  their  "universality"  that  causes 
sujjerposition. 

general  interpretation  given  at  that  time  was  that  it  was  due  to 
inertia]  range  turbulence.  This  was  because  their  spectra  were 
of  the  form  approximating  and  it  is  well  known  that 
PSD's  of  velocity  fluctuations  of  the  inertial  range  obey  the 
relation 

V(kJ  =  c,€'/'k-'/'  (2.1) 

where  e  is  the  dissipation  rate  and  Cj  is  a  universal 
dimensionless  constant.  In  Dewan  (1979)  the  inertial  range 
turbulence  interpretation  of  that  data  was  shown  to  be  ruled  out 
on  the  basis  of  scale  size.  An  alternative  interpretation  was 
therefore  suggested  that  assumed  that  the  fluctuations  were  due 
to  gravity  waves  and  that  these  waves  interacted  in  a  cascade 
fashion  which  gave  rLse  to  Eq.  (2. 1 )  in  a  manner  analogous  to 
turbulence. 

A  second  indication  of  evidence  for  a  wave  cascade  comes 
from  the  empirically  based  GM  model  itself  where  for  most 
frequencies  (sufficiently  larger  than  tOj)  the  frequency 
dependence  of  the  kinetic  energy  density  is  given  by  (using 
GM  symbols,  approximately  -  see  next  section). 
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=  (2.2) 

Tennekes  and  Lumley  (1972)  have  shown  that,  for  the  case  of 
inertial  range  turbulence 

=  (2.3) 

(O 


where  n  is  the  ratio  of  buoyancy  frequency  at  the  depth  in 
question  to  th-c  extrapolated  surface  value.  A  one  dimensional 
PSD  with  0)  or  k  dependence  is  obtained  by  integration  over 
the  entire  range  of  wave  values  of  the  alternate  variable  as  will 
be  illustrated  below.  The  atmospheric  model  uses  n  set  to 
unity. 

TTie  New  PSD  Model 


where  Vg  is  the  PSD  for  kinetic  energy.  Eq.  (2.3)  is  obtained 
by  dimensional  analysis  (c.f.  Bridgeman  (1963)). 

If  the  e  dependence  in  Eq.  (2.3)  and  (2.1)  were  to  be 
verified  for  the  case  of  internal  gravity  waves  it  would 
constitute  further  evidence  for  the  postulated  cascade. 
(Meanwhile,  the  scalings  are  suggestive).  Unfortunately  there 
are,  at  present,  no  measurements  of  to  or  k^  spectra 
simultaneously  with  e  measurements. 

It  is  well  known  that  for  turbulence  the  cascade  goes  from 
small  to  large  k  and  (o,  the  mechanism  is  vortex  stretching, 
(Tennekes  and  Lumley,  1972),  and  interactions  are  local  in  k 
and  (0.  In  the  strongly  interacting  waves  presently  under 
discussion  the  description  of  Phillips  (1966)  p.  178  seems 
appropriate.  For  steep  enough  internal  waves  he  describes 
interactions  that  can  be,  to  use  his  word,  "promiscuous",  and 
hence  possibly  local.  At  this  time  thex  is  no  proof  that  local 
cascades  in  k  and  to  are  ruled  out.  I'he  phenomenological 
nature  of  the  model  to  be  proposed  below  obviates  the  need  to 
describe  the  physical  interaction  mechanisms  or  even  the 
directions  of  the  cascades  thus  leaving  these  topics  for  future 
discussion. 


Using  the  dimensional  argument  of  Dewan  and  Good 
(1984,  1986)  and  the  extension  of  Smith  et  al  (1987)  we  see 
that  t  =  3  in  A  ([3^),  thus 

a(pJ  =  2(i+PJ'  (4.1) 

in  the  GM  formalism.  This  leads  to  the  assumption  (c.f.  the 
last  two  references). 

\(t  (k  )  =  (l  +  P  ) 
u 3  ^  (4.2) 

* 

where  a  is  a  universal  constant  of  order  unity  which  should  be 
about  1/6  as  Smith  et  al  (1987)  showed.  In  a  similar  way  the 
conclusions  of  Dewan  ( 1979)  lead  to  the  assumption 


/  2/3  /,  „ 

a  e  ( 1  +  P, ) 
w  (k  )  = - -> — — 

x) 

V 


The  GM  (72, 75)  Formalism 

A  proper  review  of  this  formalism  is  not  possible  here 
within  the  space  permitted;  therefore,  the  reader  may  need  to 
consult  Garrett  and  Munk  (1972,  1975).  The  heart  of  their 
formalism  is  E  (to,  k),  an  energy  density  sptectrum  assumed 
separable.  [Fritts  and  Chou  (1987)  gave  evidence  of  in¬ 
complete  separability,  but  since  the  frequency  spectral  slopes 
did  not  have  a  wave  number  dependence,  the  new  predictions 
of  the  present  approach  would  not  be  critically  affected.  In 
any  case,  this  matter  is  beyond  the  scope  of  this  paper.]  GM 
72,75  uses 


The  frequency  spectra  are  based  on  E  B  (co)  which,  in 
analogy  to  Eq.  (2.3),  is  given  by 

EB(w)  =  ?-^.  (4.4) 

(0 

This  is  obtained  from  a  similitude  argument  based  only  on  the 
assumption  that  the  waves  interact  so  as  to  exhibit  a  cascade  in 
(o  space.  In  other  words,  one  assumes,  as  in  the  case  of 
inertial  range  turbulence,  that  there  is  a  cascade.  Then  EB(co) 
must  depend  solely  on  e  and  to.  Eq.  (4.4)  results  immediately 
from  forcing  both  sides  to  have  the  same  dimensions.  On  the 
basis  of  the  GM  72,  75  formalism. 


where  k^  is  the  dominant  wave  number,  P  =  k/k^,  and 

A  (p)  s  (t  -  1)  (1  +  p)~*,  where  E'B(a))  and  t  are  determined 
by  experiment.  The  total  range  integrals  of  A(p)  and  B(co)  are 
defined  as  unity.  Below,  both  t  and  EB(co),  will  be  derived 
from  dimensional  analysis.  EB((o)  is  a  dimensional  energy 
spectrum  per  unit  volume  as  compared  to  the  dimensionless 
form  E'B((o)  in  GM  72,  75. 

A  crucial  part  of  their  formalism  is  the  ratio  of  vertical  or 
horizontal  velocity  fluctuations,  v  or  u  respectively,  to  total 
energy  density  E((oJc)  given  by  the  so-called  "wave  functions" 
(derivable  from  tlie  polarizoilun  relations),  defined  respectively 
by 


~2  2  V®+<ftj*n 


~2  ((o^-(a^)«n 

and  V  (to)  s  2^ - ^ - (3.2) 


(4.5) 


(4,6) 


(4.7) 


The  e  dependences  in  Eq.(4.7).  (4.6)  and  (4.3)  are  new 
predictions  for  the  gravity  waves.  Table  1  contains  a  summary 

of  predictions  of  this  approach  including  the  values  of  k^^  and 

* 

k^  .  The  derivation  of  the  latter  can  be  found  in  Dewan 
(1990). 
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Table  1.  Proposed  Model 


k^ 

% 

(m^/s^)/ 

3/3 

V  %  /  k 

(rad/m)  (m^/s^)/(r 

k 

ad/m) 

'Fv(k)  =  j( 

a  =  !tafV“ 

_ L  l  . 

1 _ 

1 

(0 

(m"/s^) 

r 

Ir  -  <2  V 

z,  z 

/  2  2'V  y  U  - 

1(0  +  CO  / !  a  £  1(0  -  (o.  /  i 

-  '  ''  2  1  .  ,2  1 

to  :  (0  1  N  1 

/(rad/s)  (m^/s-)/(''3(l/s) 

^^(oJ/E  V  ((o')/e 

and  are  universal  constants  of  order  unity  and 


a"  =  (a  sf)),  a'  =  (3/8)  a/(s;) 


-5/ 


Details  are  in  Dewan  (1990) 


In  order  to  test  the  model  (Table  1 )  simultaneous  measure¬ 
ments  of  e,  Vu  V  (to),  Vy  y  (k^),  Vm,  (k,).  and  from  them 
k  and  k  are  needed.  In  the  ocean  sonar  would  be  used 
whereas  radar  and  lidar  would  be  the  instruments  of  choice  for 
the  atmosphere. 

Concerning  Quantitative  “Reasonableness"  of  the  Model 


Due  to  the  fact  that,  as  of  this  writing,  there  does  not  exist 
a  single  simultaneous  measuremen  of  ipy  ^  (k^)  or  ^  (to) 
together  with  suitably  averaged  e,  it  is  not  possible  to  subject 
the  model  to  immediate  test.  The  only  question  that  can  be 
addressed  right  now  is  whether  or  not  the  predictions  are 
witl.'in  reasonable  bounds  of  known  measurements.  This 
section  will  study  this  question.  Note  the  citations  used  for 
"reasonable"  parametric  values. 


Atmosphere 


The  parametric  values  used  here  will  be 
C  =  5  X  10~^  m^/s^  Barat  (1975) 


N  =  10  ^  rad/s 

-4 

(0  =  10  rad/s 

t 


Chosen  to  fit  the  relevant 
altitude  and  latitude 


a  =  0.24  -  using  above  in  Table  1 

-  0.2  -  fit  to  experiment 


Occm 

e  =  10-*  (m2/s3)  Woods  (1974),  GM  (1972,  p.  251) 

N  =  3.5  X  10  3  (rad/s)  GM  (1972),  2  (cy./hr.) 

<Oj  =  7  X  10"^  (rad/s)  GM  (1972),  0.04  (cy./hr.) 
a  =  0.256  using  above  in  Table  1 

Sj  =  0.2  (from  previous  section) 

Table  2  presents  a  comparison  between  theory  and  cited 
experiiiiental  values. 

Concerning  the  choice  of  the  value  of  e  for  the  oceanic 


case  above,  it  is  the  value  cited  for  the  upper  thermocline  and 
is  therefore  appropriate  for  the  case  n  =  1  in  the  GM  model 
used  for  comparison  in  Table  2. 

Table  2. 


Quantity 

(units) 

Theory 

Atmos. 

Exper. 

Atmos. 

Tlieory 

Oceanic 

Exper. 

Oceanic 

(2  (iTl) 

22 

40(0 

0.48 

1.0(4' 

(2  n/k2^)(km) 

7 

5(2) 

107 

170(3) 

Vu  <“i) 
|(m2/s2)/(Hz)| 

1.4x106 

10-“' to 
2x106(3' 

62 

56(6) 

Vu(N) 

|(m2/s2)/(Hz)l 

69 

63  to 

3.2x10-3(3) 

1.2x10-2 

1.7x10-2(6) 

(1)  Dewan  and  Good  (1986),  (2)  Smith  et  al  (1987),  (3)  Van 
Zaiidi  (1982),  Figure  1.  (4)  Woods  (1974;,  ^3;  noiioway 
(1981)  p.  58,  (6)  GM  (1972),  Figure  4,  and  graphed  value 
times  (N/(27t)). 

Vertical  Coordinate  Dependence 

The  average  value  of  e  as  measured  in  the  atmosphere 
depends  on  altitude.  For  example  Barat  (1975)  found 
stratospheric  values  of  e  to  be  in  the  range  KT-’  to  5  x  lO”^  m* 
s”^.  We  shall  use  the  upper  limit  in  what  follows.  At 
mesospheric  heights  (86-88  km)  Hocking  (1988)  lound 
average  values  of  £  to  be  around  0. 1  m-  s'^  From  the  model 
derived  above,  (co)  is  directly  proportional  to  e  and  hence 
one  would  expect  that  it  would  be  of  order  200  times  larger  at 
86  km  than  at  1 1  km.  The  question  raised  therefore  is  whether 
or  not  such  differences  in  y^  (to)  are  in  fact  measured  in  the 
atmosphere.  According  to  Balsley  and  Carter  (1982)  "In  the 
portions  of  the  spectra,  the  spectral  energy  density/unit 
mass  is  about  250  times  smaller  in  the  troposphere  (8  km)  than 
in  the  mesosphere  (86  km)". 

From  this  one  can  conclude  that  the  new  prediction  that 
y„  (to)  will  be  proportional  to  e  as  the  latter  changes  with 
altitude  is  reasonable. 

Regarding  the  oceanic  case,  it  can  be  shown  that  the 
present  model  implies  that  e  should  scale  as  e-y^  vvheic  y  is 
depth  and  b  ~  1.3  km  where  use  is  made  of  GM  72,  75  plus 
WKB  scaling  which  implies  that  total  energy  per  unit  volume 
scales  as  N  (see  Van  Zandt  and  Fritts  (1989)  for  example).  I 
do  not  know  if  this,  in  fact,  has  ever  been  observed.  It  is  a 
prediction. 


Conclusions 

In  this  paper  a  similitude  model  has  been  suggested  which, 
if  experimentally  validated,  would  phenomenologically  extend 
the  GM  72,  75  model.  At  the  same  time,  it  would  increase 
somewhat  our  understanding  of  the  underlying  physics  of  the 
processes  in  question.  This  model  is  based  on  a  turbulent-like 
wave  field  picture.  Wave-like  properties  (e.g.  polarization 
relations,  etc.)  are  imagined  as  coexisting  with  strong  mode- 
interactions  giving  rise  to  cascade-like  transfers  in  the  k^  and 
to  scales  and  saturation-like  proces.sps  in  k^  scales. 

A  process  known  as  "buoyancy  range  turbulence"  has 
been  described  in  the  literature,  sometimes  in  cormection  with 
GM  72,  75.  See  Luniley  (1964),  and  Wtinstock  (1985).  It  is 
possible  that  the  model  offered  h^.e  will  turn  out  to  be  an 
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improved  description  of  such  "buoyancy  range  turbulence", 
but  since  wave-like  properties  are  assumed,  the  word 
turbulence  would  have  to  be  modified.  "Wavelike  turbulence" 
might  be  appropriate  (compare  Gage  (1979)). 

Finally  it  should  be  mentioned  that  phenomenological 
models  are  helpful,  not  only  to  applied  scientists  and 
engineers,  but  also  to  theoreticians.  It  is  therefore  hoped  that 
this  work  will  be  of  assistance  in  the  development  of  a  "first 
principles"  physical  theory  of  GM  72, 75. 
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